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We study finite- frequency quantum noise and photon-assisted electron transport through a wide 
and ballistic graphene sheet sandwiched between two metallic leads. The elementary excitations 
allow as to examine the differences between effects related to Fabry-Perot like interferences and 
signatures caused by correlations of coherently scattered particles in electron- and hole-like parts of 
the Dirac spectrum. We identify different features in the current- current auto- and cross-correlation 
spectra and trace them back to the interference patterns of the product of transmission- and reflec- 
tion amplitudes which define the integrands of the involved correlators. At positive frequencies the 
correlator of the auto-terminal noise spectrum with final- and initial state associated to the measure- 
ment terminal is dominant. Phase jumps occur within the interference patterns of corresponding 
integrands, which also reveal the intrinsic energy scale of the two- terminal graphene setup. The 
excess noise spectra, as well as the cross-correlation ones, show large fluctuations between positive 
and negative values. Oscillatory signatures of the cross-correlation noise are due to an alternating 
behavior of the integrands. 
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INTRODUCTION 



Ballistic electron transport ^'^ in two-terminal 
graphene systems is in the focus of intensive stud- 
ies ever since the pioneering experiments on single-layer 
carbon. ^'^ The Dirac Hamiltonian^'^ describes charge 
transport close to the charge-neutrality point and 
leads to a linear dispersion relation e = hvph. This 
allows to observe several relativistic phenomena in 
solid-state system, such as Klein tunneling^^^ or 
the Zitterbewegung.^^^^ In the very early works on 
graphene the minimal conductivity^^'^^'^^'^^'^^ G ^ e^ /h 
per valley and pseudo-spin at the charge- neutrality 
point has been found and stimulated the research on 
current and noise properties. The current-current 
correlations around the minimal conductivity lead to 
a zero frequency sub-Poissonion Fano factor with a 
maximal value of F = 1/3,^^^^ remarkably similar to 
diffusive systems as disordered metals. ^^^'' The sup- 
pression of the Fano factor below the Poissonian value 
originates from noiseless, open quantum channels that 
are found at all conductance minima in graphene-based 
two-terminal structures'^, and can be explained as an 
interplay between Klein tunneling, resonant tunneling 
and pseudo-spin matching. This pseudo-diffusive behav- 
ior'^ is due to the special band-structure of graphene. 
Without impurity scattering, coherent transport through 
such a graphene sheet'^ gives rise to the same shot noise 
as in classical diffusive systems. The opening of a gap^^ 
in the quasiparticle spectrum leads to an enhanced 
Fano factor.'^ Such a gap can be opened for example 
by photon-assisted tunneling, as shown recently for the 
case of a graphene p-n junction^' with a linear potential 
drop across the interface. ^^'^^ There, Landau-Zener like 
transitions stimulated by photon emission or absorption 
via resonant interaction of propagating quasiparticles 
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FIG. 1: Wide graphene strip (W ^ L) sandwiched be- 
tween two heavily-doped, metallic graphene leads. The Fermi- 
level of the sheet can be tuned by a center gate voltage Vg. 
Electron- and hole states are injected via the time-dependent 
bias- voltages in left- and right leads /xl/hI^)- 



in graphene with an irradiating electric field lead to 
hopping between different trajectories. 

The scattering approach as put forward by Landauer 
and Biittiker^^ has been applied to ac-driven charge 
transport^^^^ through a metal- graphene interface with 
an abrupt potential change. ^^ The metal can be formed 
by a graphene lead strongly electrostatically doped 
by a gate potential, thus shifting the Dirac point far 
away from the Fermi energy. In this work we adopt 
the formalism and parameterization introduced in 
Refs.^^'^' and calculate the finite-frequency current- 
current correlations at finite dc- and ac-bias voltages in 
the system depicted in Fig. 1. We complement recent 



results on ac-transport in Fabry-Perot graphene devices 
of Ref.^^, in which the influence of different boundary 
conditions, i.e. zigzag or armchair configurations, on the 
Fabry-Perot patterns in a combined Tien- Gordon/tight- 
binding approach has been investigated. The influence 
on transmission properties of a time-dependent potential 
barrier in a graphene monolayer has been investigated.^^ 
In our work the transverse boundary effects are central 
and we assume so-called infinite mass boundaries^^'^^ 
describing a short but wide (L <C W) graphene strip. 
We focus on the interplay between the Dirac-spectrum 
with the Fabry-Perot interferences. 

Interestingly, the well-known oscillations as function of 
gate voltage on a scale of the return frequency hvp/L^ re- 
lated to the length L of the graphene sheet, can be seen as 
a reminiscence of Zitterbewegung.^^ The role of the com- 
plex reflection amplitude and the onset of contributions 
of scattering states coming from terminal a and being 
scattered into terminal /3 will be the key characteristics 
in our discussion of the results for the noise as function 
of bias voltage and frequency. As a consequence of these 
onsets the oscillations add up de- or constructively de- 
pending on the precise values of voltage and frequency. 
In our setup, the separation of oscillations caused by the 
Fabry-Perot reflections and effects caused by the band- 
structure of the Dirac Hamiltonian is a priori not obvious. 
In both cases phase-coherent transport is essential. How- 
ever, for charge injection either into the conduction or the 
valence band only, effects like Zitterbewegung should not 
be present and all oscillating features of the noise spectra 
have to be of Fabry-Perot nature. 



II. DIRAC EQUATION AND SCATTERING 
FORMALISM 
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FIG. 2: (color online) Transmission probability T(e) = |t(e)| 
as a function of energy and transverse momentum q. 



Making use of the Tien-Gordon ansatz, we write the so- 
lution to the time-dependent Dirac equation as a sum 
over PAT modes: 
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The advantage of this ansatz is that the scattering prob- 
lem has to be solved for the time-independent case only. 
Therefor, in terminals 7 = L, i? we define stationary so- 
lutions ^o(x, t) = ^e(x)e~*^^ by the equation 



The ballistic graphene^^^^ sheet considered in the fol- 
lowing can be described by the two-dimensional Dirac 
equation for the two-component spinor ^ = (^1,^2)^ 
with indices referring to the two pseudo-spins of the car- 
bon sub-lattices. Throughout this work we will neglect 
inter- valley scattering and Coulomb interactions. We 
only consider the interaction of the electrons with the 
radiation field in the form of photon-assisted transitions. 
With Fermi velocity vp the Dirac equation can be cast 
into the form 
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The electrochemical potential /i(x, t) includes static and 
harmonically driven potentials in the leads plus a static 
gate voltage in the graphene sheet. 
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The basis states in graphene can be constructed as a 
superposition of left- and right movers. 
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a{e) describes the angle between the momentum of a 
quasiparticle and it's y-component q in region x = . . . L 
of the graphene sheet. Then the pseudo-spinors can be 
parametrized as 
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Here the dispersion is given by e = hvF\^q'^ + k'^- 
wave vector k{e) and the angle a{e) are defined as 



a{e) = arcsin 
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Therewith, and neglecting transverse momentum due to 
high doping, we have the basis states 
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in the leads. Additionally shifting the Fermi surface of 
the graphene sheet away from the Dirac point, and thus 
changing the concentration of carriers, is incorporated 
into the gate voltage eVg. For |e + eVg\ < \hvFq\ we have 
evanescent modes, ^^ with imaginary a{e) and k{e). Oth- 
erwise we have propagating modes and scattering is only 
at X = 0, L. 

Irradiating the two-terminal structure with a laser^^ can 
be described by a harmonic ac-bias voltage with driving 
strength a = eVac/^ as discovered in the pioneering 
paper by Tien and Gordon^^. Their theory can be in- 
corporated into the scattering formalism^ ^'^^ and we are 
applying it here to the two-terminal graphene structure. 
We take the two valleys and two pseudo-spin states of 
the carbon lattice into account in the pre-factor of the 
current operator of reservoir 77, which reads 
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dede / dqJi {a^) J^ {as) 



X d\{e - lhw)A^s{r], e, e)ds{e - khw)e'^'-''^*^^ . (15) 

Indices 7, S run over reservoirs L, R. Summation over all 
modes of y-momentum is replaced by an integral since 
W ^ L. Scattering is contained within the current ma- 
trix A^§{r]^ e, e') = 5j^^5j^s — s^ {e)Sj^s{e') of a current be- 
tween leads 7 and 5 measured in lead rj via the energy- 
dependent scattering-matrix 
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t{e) r'{e) 



(16) 



The scattering matrix connects in- and outgoing scatter- 
ing states at the two barriers and is calculated in Ap- 
pendix B by matching the wave functions at x = 0, L. 
Here we write the results for refiection and transmission 
amplitudes in an alternative version: 
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FIG. 3: (color online) Left: Conductivity a{uj,a) — 
{L/W)G{uj^a) as a function of dc-voltage applied across the 
two-terminal setup. We show curves for various ac-driving 
strengths a applied to a) both reservoirs (a = 0) and b) to 
the left reservoir only (a = 1) with ujL/vf = 5. 



We assume identical scattering for quasi-particles inci- 
dent from left and right, so t{e) = t' {e) and r(e) = r'{e). 
r(e) vanishes if k{e) = 7rn/L, with integer n. The corre- 
sponding modes in ^/-direction are determined by 
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giving rise to special features of the current fiuctuations, 
going along with the phase jumps of ttL/Hvf in r(e) we 
discuss later on. At the Dirac point transmitted quasi- 
particles at perpendicular incidence perform Klein tun- 
neling via evanescent modes, leading to finite transmis- 
sion probability T(e) = t^{e)t{e) at small transverse mo- 
mentum, see Fig. 2. 



III. DIFFERENTIAL CONDUCTANCE 

Since the average current has only a zero-frequency 
component, PAT events in the conductance^^'^^ can only 
be studied by inducing photon-exchange via a time- 
dependent voltage as it is, for example, generated by 
irradiating the setup with a laser beam. Different po- 
larizations of the coupled light field lead to different ac- 
driving in left and right leads. Such an asymmetry can 
be described by a parameter a G [—1,1] which varies the 



driving in the leads via aL/R 
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FIG. 4: (color online) Schematic view of the different regions 
which occur in the integrands of the correlators contributing 
to the finite- frequency shot noise spectrum. 



hidden behind the contribution from the terminal where 
driving gets small, as in Fig. 3 b) with a = 1. From the 
oscillations with period proportional to L, we expect no 
measurable effect on conductivity or shot-noise^^'^^, as 
in the case without ac-driving and for the zero-frequency 
Fano factor. In the scattering approach they are simply 
because the transmission function oscillates as a function 
of energy. But imperfections of real samples, as impuri- 
ties^"^ or lattice-mismatch, lead to scattering events. Due 
to this randomizing effect on the path-lengths for prop- 
agating quasi-particles the calculated oscillations are av- 
eraged out in experiment^^'^^. 



IV. FREQUENCY-DEPENDENT SHOT NOISE 



call the driving symmetric (in the amplitudes Kic,l/r) if 
a = and asymmetric if a = ±1. For clearness we will 
only discuss a = 0, ±1 since intermediate values are just 
a mixture of those limiting cases. For arbitrary a the dif- 
ferential conductance can be derived from Eqn.(15) by 
taking the statistical- and time average and differentiat- 
ing with respect to voltage. At /c^T = it reads 
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Different orders m of PAT do not mix but have to 
be summed up resulting in independent contributions 
G^(cc;,Qf) to differential conductance. Since G{uj^a) only 
depends on the Besselfunctions squared, these pre-factors 
will always be positive. The influence of the driving 
strength a on conductivity (j{(jJ^ a) = {L/W)G{(jJ^ a) as a 
function of dc-bias is plotted in Fig. 3. PAT events lead 
to a substantial enhancement of the conductivity around 
zero dc-bias, because more channels are available in com- 
parison to the case without time-dependent voltages. At 
large dc-bias voltages this effect gets negligible since the 
transmission probability of the graphene sheet, see Fig. 2, 
is not vanishing at large energies. Thus, those contribu- 
tions built a dominant background. Conductance at ar- 
bitrary dc- and ac-bias is a sum of two integrated trans- 
mission probabilities, where the integrand exhibits cross- 
ings of the two independent interference patterns, as in 
region IIB in Fig. 4 a). Each Gmip^) shows a transi- 
tion from a region with an oscillating, but in average not 
increasing contribution to conductance for dc-bias volt- 
ages \eV/2\ < \mhuj\, to a regime with a linear increas- 
ing background at larger dc-bias voltages. The photon- 
energy mhw introduces a phase shift in the oscillations of 
Gm{^-, (y) as a function of dc-bias voltage, so for different 
m we can have local minima or maxima at eV = 0. After 
summation, conductivity can also show a local minimum 
or maximum at eV = 0, as it can be observed for the vari- 
ous values of a in Fig. 3 a). If |a| tends to one this effect is 



To get full informations on current-current correla- 
tions we study the non-symmetrized noise-spectrum as it 
can be detected by an appropriate measurement device 
in the quantum regime. ^^^^ 

We allow harmonic ac-driving eVac cos{ujt) in the leads, 
so in Fourier space the current-current correlations are 
defined as 
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The non-symmetrized shot noise correlates currents at 
two times: 



Sap{t,t\w) = {AL{t)AIp{t') 



(22) 



with variance AIa{t) = Ia{t) — {Ia{t))- Of experi- 
mental interest are the fluctuations on timescales large 
compared to the one defined by the driving frequency 
uj. Thus, as in^^, we introduce Wigner coordinates 
t = T -\- r/2 and t' = T — t/2 and average over a driv- 
ing period 2i\ juj. Then, the noise spectrum is defined by 
the quantum statistical expectation value of the Fourier- 
transformed current-operator /^(r^) via Sa,^(Vt^Q! ^uj) — 

2i^Scc^(Sl,uj)^{Sl ^ ^) = (/a(^)//3(^0)- Safsi^^oo) is 
nothing but the Fourier transform of Sai3 (t, cj) . Similarly, 
in the case without ac-driving the noise is only a func- 
tion of relative times r = t — t' . In order to keep notation 
short, in the dc-limit we write SajsiS^) '•= Sai3{^^oj = 0). 
To get a deeper insight into the underlying processes of 
charge-transfer we split the noise into four possible cor- 
relators^^, defined by 



Sll{^,^) '= ^ Cc,^/3(^,CJ). 
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a,(3=L,R 



The correlators itself can be seen as the building- 
blocks of noise spectra where different combinations de- 
scribe noise detected by corresponding measurement se- 
tups. ^'^'^^ First we discuss Sll{^) := Sll{^,^ = 0), the 
case when no ac-driving is present. We also skip uj in the 



arguments of the correlators. Then evaluation of Eqn. 
(21) at ksT = leads to the expressions: 



transmission- or reflection probabilities: 
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At flnite dc-bias voltages correlations with initial and fl- 
nal state related to the measurement terminal L are spe- 
cial in the sense that they can not be written in terms of 
probabilities at flnite frequency. For symmetrized noise, 
Biittiker^^ discussed the essential role of the complex re- 
flection amplitudes in elastic electron transport and how 
they determine the equilibrium current fluctuations. In 
the quantum regime at /c^T = 0, the equilibrium fluctu- 
ations are given by 
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For flnite dc-bias the reflection amplitudes entering 
Cl^l(^) play the same essential role as in equilibrium, 
in the sense that flnite-frequency current fluctuations are 
non-zero even for vanishing transmission. The combina- 
tion of scattering-matrices of the correlators integrands 
which enter in the current-current cross-correlation spec- 
trum 
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are substantially different than in the ones for the auto- 
terminal noise. Most of all, at flnite frequency none of 
the complex correlators can be written as an integral over 
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Unlike for symmetrized noise, quantum noise^^'^^ spec- 
tra discriminate between photon absorption (Q. > 0) and 
emission {Q < 0) processes between quasi-particles in 
graphene and a coupled electric fleld^'''^^'''^ ''^. Energy 
for photon emission has to be provided by the voltage 
source, so at ksT = the Heaviside-Theta functions en- 
sure that only terms satisfying this condition contribute 
at negative frequencies. In the dc-limit, our choice of 
chemical potentials —jul = /^i? = eV/2 > and the 
fact that the measurement is performed at reservoir L, 
leaves only C^^i;,(^) ^ if (^ < 0. When additional ac- 
voltages are present none of the correlators of Eqn.(Al) 
is given in terms of probabilities and integration bound- 
aries are changed by :^mhw. Then all correlators can 
contribute at frequencies 1] < 0. 



V. QUALITATIVE DISCUSSION 

A good starting point to interpret results for conduc- 
tivity and shot-noise spectra is to examine the involved 
integrands in Eqs. (24) and (27). Figure 4 provides a 
schematic overview of the different regions occurring in 
the 2D-plots of Figs. 5-10. We show the real parts of 
integrands either as a function of {q^ e) as in scheme 4 
a) or of ((7,1^) as in scheme 4 b). The former is divided 
by the four envelopes q = \e\ and hvFq = \e -\- hQ\ into 
six areas: I, where the regimes 11^ and II5 of evanescent 
modes are merging and the areas 111^,111^, III^ of propa- 
gating modes. Area III5 is deflned by the two lines with 
origins (g' = 0,e = 0), (g = 0, e = —HQ) and intersec- 
tion {hvpq = h£t/2^e = —hQ/2). Areas in scheme 4b) 
are separated by hvFq = |e + hQ\ and the dashed hori- 
zontal line hvpq = \e\. The transmission probability flts 
into this scheme when the horizontal separation is absent 
so we are left with areas 1^ and 2^1/5. Then area 1^ in- 
cludes the black region of Fig. 2 where no transmission is 
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FIG. 5: (color online) Real parts of integrands of the 
four correlators (Eqn. 24) contributing to the shot-noise, 
namely a) 0.25|1 - r'^ieyie + ;^Q)|^ b) T{e)T{e + HQ), c) 
R{e)T{e + hn) and d) R{e + hQ)T{e). Here the energy is fixed 
e = corresponding to vanishing dc-bias. The correlator 
in a) cannot be written in terms of a probabilities, except 
in the zero frequency limit the integrand results in T^(e). 
Correlator b) contains one transmission probability at zero 
energy that is only non-zero at small q. Since for small 
transversal momentum R{e) decays as q~^ the correlator 
c) tends to zero in this regime and otherwise mimics the 
behavior of T(e). Integrand d) is also restricted to low 
transverse momentum because T(e) = otherwise. 



possible, and the regime of evanescent modes with finite 
transmission probability for small |e| < hvF\q\ around 
e = due to Klein tunneling. In regimes 2a/b a hyper- 
bolic shaped interference pattern with oscillations along 
e is prominent, where the period of oscillations is on the 
order of hvf/L for small HvfIqI <^ |e|- Figure 5 shows the 
relevant integrands of the four correlators Ca^f3{^) con- 
tributing to the finite frequency quantum noise, plotted 
as a function of {q^ Q) when e = 0. Then the imaginary 
part of r*{e)r{e + hQ) leads to finite contributions in the 
region I^ and I5 in figure 5a). T(e = 0) is only non-zero 
for small g, so integrands b) and d) vanish for large q. 
Since R{e) = 1 — T(e), integrand c) vanishes when q ^ 
and otherwise resembles the shape of T(e). 
Finite e, as in Fig. 6, introduces another interference pat- 
tern for propagating modes. In region 1^ non-zero val- 
ues are possible and in 2^ and 25 the usual interferences 
occur. For g'-values below hvF\q\ = |e| this additional 
pattern can be seen in region I5. The interplay of both 
patterns leads to phase jumps of nL/hvp in regions 3a 
and 85. These phase jumps can be determined by re- 
quiring |r*(e)r(e + hVt) -1^ = 1 in Eq.(24a), Fig. 6a). 
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FIG. 6: (color online) Real parts of integrands of correlators 
(Eqn. 24, see also Fig. 5) contributing to the shot-noise for 
fixed energy eL/hvp = 20. At finite e there is an additional 
interference pattern along q if hvF\q\ < |e|, leading to phase 
jumps in the Integrand of correlator a), the one where initial 
and final state belong to the measurement terminal L. When 
the integrand can be written as a product of probabilities, 
see b)-d), the phase jumps are absent but two independent 
interference patterns are found. 
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FIG. 7: (color online) Real parts of integrands of the correla- 
tors (Eqn. (24), see also Fig. 5) contributing to the shot-noise 
with fixed frequency QL/vf = 20. Analogous to Fig. 6 but 
as a function of (q, e). Phase jumps occur in the intervall 
—hQ < e < in integrand a), region 111^ of Fig. 4a). The in- 
terplay of the two interference patterns can also be observed 
at larger energies and transverse momenta for hvF\q\ < |e|, 
hvF\q\ < |e + hQ\ in all integrands a)-d). 




-40 -20 20 40 -40 -20 20 40 




-40 -20 20 40 




0.5 



20 40 




-40 -20 20 40 

e[fiYp/L] 



-40 -20 20 40 

e[ftyp/L] 




0.5 



e[fiYF/L] 



-40 -20 20 40 

e[fiYF/U 



-40 -20 20 40 

e[fiYF/U 



FIG. 8: (color online) Real parts of integrands which ap- 
pear in Eqn. (27) contributing to the cross-correlation shot- 
noise for fixed energy eL/hvp = 20 (top) and fixed frequency 
QL/vF = 20 (bottom), namely a),c) 3^[t*(e + ;iQ)t(e)(r*(e + 
hQ)r{e)-l)] andb),d) 43f^[r*(e)t(e + ;iQ)r*(e + ;iQ)r(e)]. Due 
to symmetry reasons the integrands are identical when in- 
terchanging index labels L^R. As a function of frequency 
integrand a) leads to strongly oscillating contributions to the 
noise spectrum. These oscillations are reduced due to the 
alternating behavior along q in cross-terminal contributions 
b). In c), d) the integrands are plotted as a function of (q, e) 
where they reveal a similar structural difference. 



FIG. 9: (color online) Real parts of integrands appearing in 
the correlators of Eqn. (Al) contributing to the shot-noise for 
driving frequency ujL/vf = 7.5 and fixed frequency QL/vf = 
20. The integrands are a) 0.25(1 - r*(e)r(e + /iQ))(l - r*(e + 
hQ + nw)r{e + nw)), b) t*(e)t(e + nQ)t*(e + ;iQ + ;icj)t(e + ;kj), 
c) r''{e)t{e + /iQ)t*(e + /iQ + /kj)r(e + huu) and d) ^"{€^{6 + 
hQ)r'' (e-\-hQ-\-huj)t{e-\-huj). When two frequencies are present 
none of the correlators can be written in terms of probabilities 
and additional phase jumps come into play. 



Therefor r^{e)r{e + hQ) has to vanish, what is fulfilled 
by the transversal momenta of Eq. (19). The condition 
|r*(e)r(6 + hQ) — Ip = 4 for a maximum in the inte- 
grand leads to modes which experience Klein tunneling. 
Actually, this correlator can be written as integral over 
1 + R{e)R{e + HQ) - 2[R{e)R{e + HQ)]^^'^ cos(<l>(e, n)) in- 
cluding a scattering-phase ^(e, ft) = Arg [r*(e)r(e + Q)]. 
Thus it describes events containing the scattering-phase 
between time-reversed paths of electron-hole pairs sep- 
arated by the photon energy hft reflected back into the 
measurement terminal. The effect of the phase shifts on 
the integrands interference patterns is also obvious in the 
(g', e)-plot of Fig. 7 a), region III5. Figs. 7 b)-d) show a 
similar interference pattern although the corresponding 
correlators are defined in terms of probabilities. 
Concerning cross-correlation noise, the integrands occur- 
ring in Eqs.(27) show alternating patterns of positive and 
negative values. The ones which describe auto-terminal 
contributions to Slr{^) (Eqs.(27a) and (27b)), as in 
Fig. 8a), have an alternating sign along Q. In the cross- 
terminal ones (Eqs.(27c) and (27d)), as in Fig. 8b), 
the additional interference pattern along q introduces an- 
other change of sign. Plots 8 c) and d) show a similar 
behavior as functions of (q^e). When ac-bias voltages 
introduce the driving frequency cj, the integrands struc- 
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FIG. 10: (color online) Real parts of integrands of the cor- 
relators contributing to the shot-noise with driving frequency 
ujL/vf — 7.5 and fixed energy eL/hvF — 20. Analogous 
to Fig. (9), but as a function of {q^^). As a consequence 
of PAT horizontal interference lines occur for transverse mo- 
menta hvF\q\ < |e + 7nhuj\ as in figures a), c). 



tures become even richer but also less clear, as in Fig. 9 
and Fig. 10. Then alternating signs in all contributions 
to auto-correlation noise are observed, except for the cor- 
relator with initial and final sates in the measurement 
terminal. This results in peculiar oscillatory features in 
the interference patterns at combinations of all involved 
energies e, hQ^ mhuj. Predicting the effect of such features 
on the noise spectra from the plotted integrands is then 
almost impossible because one still has to average over 
all possible energies and g-values by integration. 



VI. AUTO-CORRELATION NOISE 

In contrast to conductivity, the shot-noise spectrum in 
general couples different orders of PAT events, expressed 
by the product of four Besselfunctions of arbitrary order. 
But since the driving is fixed, non-vanishing contribu- 
tions exist only up to a certain order depending on the 
precise value of a. When time-dependent voltages are 
present, current fluctuations of Eq. (Al) contain prod- 
ucts of four scattering matrices, each with a different en- 
ergy argument. After performing the dc-bias limit only 
transitions between e and e + hQ are left. 



A. Shot noise spectrum 

In the regime eV^ hVt^hw <C hvp/L^ the scattering 
matrix can be treated as energy-independent. Then, 
as for a single level quantum dot in the broad-band 
limit, asymmetric quantum noise as function of fre- 
quency is the sum of four straight lines, with kinks 
at h£t = 0, ieF.^^'^^ For vanishing dc-bias we have 
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as long as 1] <C vp/L. The richer regime, when 
eF, hVt^hw > hvp/L, additionally exhibits strongly 
oscillating integrands. Those oscillations are purely due 
to propagating modes as it is also clear from interfer- 
ence patterns of the integrands in Figs. 5-9, regions 
lla,b and Ilia, 5, c- In the special case of perpendicular 
incidence {q^ a{e) = 0) we have Klein tunneling, thus 
the frequency-dependence of the correlators is linear for 
this mode. Then Ca^(3{^) = if a ^ /3 since R{e) = 0. 
Otherwise the C^^^i^) mirror the interference patterns 
of the integrands. So the noise spectrum (Fig. 11 solid, 
thick curve) shows oscillations on the scale of L/hvp 
in the regime eV^hft^hu :^ hvp/L, similar to the shot- 
noise at zero-frequency as a function of gate voltage^^. 
Although present in all four correlators, the oscillations 
show up in the noise spectrum mainly via C^^2.(^) ^^ 
the terminal where the fluctuating currents are probed. 
That is because the correlator itself as well as the 
amplitude of the oscillations are significantly larger 
than for other contributions. Therefore, in comparison 
to the absorption-branch (positive frequencies) the 
emission-branch of the spectrum (negative frequencies) 
shows only small shot-noise. Indeed all correlators 



spectrum in units of 27vh/e . We compare a setup where dc- 
bias voltages are fixed symmetrically around the Dirac point 
(top,) with the case when eVoL/hvp = 2eV/hvF = 10 (bot- 
tom). Thick lines: Shot-noise and correlators. Thin lines: 
Derivatives with respect to frequency. Contributions from 
Cl^l{^) are dominant at positive frequencies. Top: Special 
features in the derivatives are seen for frequencies HQ < eV in 
the R ^ R contribution, when the lower bound of the energy- 
integration interval approaches the Dirac point (compare to 
Figs. 6, 7). Bottom: The distance to the Dirac point is in- 
creased by the offset voltage. Therefore oscillatory features 
appear in a larger frequency interval and in all four corre- 
lators, since integration boundaries in all contributions are 
crossing the Dirac-point with increasing Q. 



except Cr^l{^) vanish when 1^ < since the energy 
for the emission of a photon has to be provided by the 
voltage source. Especially the contribution dominant 
at positive frequencies vanishes: Cj^^j^{Q) = if hQ < 0. 

We are considering the limit ksT = where the cor- 
relators integration windows are exactly determined by 
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FIG. 12: (color online) Real parts of auto-correlation excess- 
noise spectrum in units of 27Th/e^ (thick lines, upper panels) 
and derivatives (thin lines, lower panels) with dc-bias sym- 
metrically applied around the Dirac point (top) and for finite 
eVb = 2ey (bottom) . By subtracting the noise at zero dc-bias 
the divergent background is removed. The structure and es- 
pecially the oscillatory behavior are coined by auto-terminal 
contributions of Eq.(27a) related to the measurement terminal 
L. The jump in the derivative of Cl^r{^) is present because 
this correlator does not contribute for frequencies HQ < eV. 
By applying an offset voltage eVoL/hvp = 10 a complicated 
structure emerges, best visible in the derivatives. 



the chemical potentials. At finite temperature this so- 
defined onsets of the four contributions as a function 
of frequency are smeared out by the broadening of the 
Fermi-functions. Clearly a gate voltage does not affect 
these onsets since it does not enter in the Fermi func- 
tions of the leads, but it still changes the transmission 
function resulting in a modified spectrum. Those limits 
of energy-integration, as well as their position relative to 
region III5, result in features in the noise spectra besides 
the discussed oscillations. In order to clarify the role 
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FIG. 13: (color online) Real parts of auto-correlation current- 
current fluctuations in units of 2'Kh/e^ as a function of dc-bias 
for fixed frequency hTl (thick lines, upper panels). We com- 
pare the symmetric setup without dc-bias offset (top) and 
when eVoL/hvp = 10 (bottom). Thin lines (lower panels) 
are used for the derivatives with respect to voltage. Top: 
Due to symmetrically applied bias voltage the noise and the 
auto-terminal contributions are symmetric in the voltage de- 
pendence and Sa^f3{^)\Y = S/3^a{0.)\_y if a / /3. Bottom: 
By applying an offset voltage we are breaking the setups sym- 
metry. Auto-terminal terms are then symmetric with respect 
to eV = =b2eVb while the summed up noise is asymmetric. 



of the Dirac Hamiltonian in comparison to the role of 
pure Fabry-Perot interferences, we compare results when 
the charge injection is only in the conduction or valence 
band by shifting the dc-bias voltages above the Fermi 
energy of the graphene sheet via the offset voltage Vq 
in /J^L/R = ^eV/2 + cVq. C]^^j^{^) can never see the 
regime —hQ < e < when eVo = 0, as in the upper plot 
of Fig. 11. Thus the oscillations visible in the derivate 
have a well defined period over the whole spectrum on 
top of a linearly increasing background. When an offset 
voltage eVo = ^eV is applied, as done when calculat- 
ing the spectra for the lower plot of Fig. 11, Cr^r{CI) 
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shows a complicated frequency dependence for small Q. 
Contribution Cl^l(^) describes correlations of scatter- 
ing states emanating from the left reservoir reflected back 
into the same reservoir. We will discuss this contribution 
now in detail: Special features for small frequency are 
due to the interplay of the integration boundaries with 
the various regions in Fig. 4a) occurring in the integrands 
(g', e)-dependence of Fig. 7a). Integration is over all q- 
modes and from e = —eV/2-\-eVo — hQ to e = eVo — eV/2. 
When eVb = 0, eF = this corresponds to —hQ < e < 0^ 
regions III5 and partly II^^^ of Fig. 4a). Now at finite 
eV, eVo as in Fig. 11, the integration window can include 
region III5 completely, partly, or not at all, resulting in 
variations of the spectrum. At small hQ, features in the 
integrands interference patterns have stronger impact. 
This can be seen from strongly non-harmonic features 
of the noise spectrum, e.g. in Cl^l{^) and Cr^l{^) 
for eVo = 2eV. For large frequencies averaging leads 
to nearly harmonic oscillations on top of the increasing 
background. With the chosen parameters the distance of 
the chemical potential /i^ to the charge-neutrality point 
is given by e{-V/2 + Vo)L/{hvF) = 7.5. Around the 
corresponding frequency the oscillatory behavior of the 
spectrum is modified and flattened due to a reduced frac- 
tion of propagating modes. Raising the frequency further 
increases this fraction again and oscillations are roughly 
harmonic with period ttL/Jivf^ best visible in the deriva- 
tives dCL^L(S^)/dQ. of Fig. 11. That is also the point 
where the lower bound of energy integration starts to in- 
clude the special interference pattern of the integrands 
around the energy interval —Vt < e < 0, region III5. 
Cr^r{^) is not influenced by the measurement terminal 
itself, but probes transmission probabilities via scattering 
events which are related to the right terminal only. An 
analogous behavior of the spectrum as before is found, 
this time with a distance e{V/2 -\- VQ)L/{hvF) = 12.5 of 
the lower integration boundary to the charge neutrality 
point when hQ = 0. Now increasing frequency is going 
along with a decreasing slope of the derivative with re- 
spect to frequency until the Dirac point is reached. There 
the slope increases again since more open channels be- 
come available. The same interpretation also explains 
features in the interval hft < eV of the auto-terminal 
correlators shown in Fig. 11, when Vb = 0. E.g. the 
spectrum of the correlator Eq.(24b), with initial and fi- 
nal state in the right lead, exhibits a reducing slope until 
hQ = eV/2 from where on the oscillations have a well 
defined period. The dCR^R(p)/dVt curve has a maximal 
slope at h£t = eV when positive and negative energies 
with same magnitude are present. For higher frequencies 
oscillations have again a well-defined phase. 
We also study the excess noise at finite frequencies: 
5'exc(^,^) •= S{Q,uj)\^y - S{ft,uj)\^y^Q. Subtracting 
the noise at zero bias-voltage removes the divergent con- 
tributions from the noise spectrum. Then oscillating fea- 
tures due to bias-voltages are more obvious since they 
are now also prominent in the noise spectra of Fig. 12, 
not only in derivatives. When eVb = the excess noise 



(thick, black, solid curve) is purely positive for hft <C eV 
while for hQ > eV it is oscillating around zero, because 
then cross-terminal contributions Ca^isi^) cancel each 
other up to a constant offset acquired at small Q. This 
offset is compensated by the L ^ L contribution. Os- 
cillations of this contribution have again a considerable 
impact on the excess noise spectrum. In the lower plot 
of Fig. 12 the offset voltage is fixed to eVb = 2eV. For 
low frequencies hft < eV, complicated oscillations oc- 
cur in all contributions to excess noise and are accom- 
panied by a strongly increasing slope up to frequencies 
hQ > eVo-\-eV/2. As for the noise itself, the frequency of 
the oscillations is determined by hQz = 2eV and equals 
the frequency expected from the Zitterbewegung of rela- 
tivistic Dirac fermions^^. This frequency corresponds to 
a period of T = tt in our plots. It would be interesting 
to test experimentally if those much more pronounced 
oscillation, compared to the overall shot-noise, can be 
detected in spite of randomization effects of imperfec- 
tions on the quasi-particles path lengths. In summary, i) 
the impact of the Dirac Hamiltonian on the frequency- 
dependence of auto-terminal current fluctuations leads to 
peculiar oscillation for energies in the vicinity of the Dirac 
point as an interplay of Klein tunneling, phase-jumps in 
the correlators and their energy-integration limits. And 
ii) oscillations due to the FP setup have a constant phase 
for high energies when propagating modes are dominant. 
Then dS^^^ (ft) / dft oscillates between positive and nega- 
tive values with a period as it is expected from the effect 
of Zitterbewegung. 



B. Dc-bias dependence at finite frequency 



Analogous to the spectrum, the dc-bias dependence for 
fixed frequency is featureless in the regime eV^ hVt^hw <C 
hvp /L^ except the pronounced onsets of the four corre- 
lators. This is not surprising when looking at the deriva- 
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tives with respect to voltage: 
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FIG. 14: (color online) Real parts of cross-correlation spec- 
trum in units of 27rh/e^ when eVb = (top) and eVoL/hvF = 
3eVL/hvF = 18 (bottom). Without an offset voltage (eVb = 
0) auto-terminal contributions are identical, as well as cross- 
terminal ones at large frequencies hO, ^ eV. At finite eVb 
the asymmetric bias voltage is reflected in the frequency- 
dependence of the auto-terminal correlators by their different 
heights and the shift of the oscillations maxima. 



dCL^L e^ein) 



dq 



dV 'iTTh 

- |1 - r*{-eV/2 - hn)r{-eV/2)\ 

CO 



|1 - r%-eV/2)r{hn - eV/2)\^ 

(28a) 



/ 



dq 



dV Anh 

[T{eV/2)T{eV/2 + hQ) - T{eV/2 - hQ)T{eV/2)] 



(28b) 



dCL^R e^Q(Q.-eV) 



dq 



dV Airh 

[T{hn - eV/2)R{-eV/2) - T{eV/2)R{eV/2 - hft)] 



(28c) 



dV 27rh 



oo 

/ 



dq 



[T{eV/2)R{eV/2 + HQ) - T{-eV/2 - hn)R{-eV/2)] 

(28d) 

Scattering amplitudes are roughly constant for a given q- 
mode in this regime, then correlators are straight lines as 
afunctionofdc-bias voltage. E.g. a special situation that 
could exhibit interesting physics is when some derivatives 
are zero. But this is, due to symmetry arguments, only 
possible at eV = 0, ±hQ^ proofing a zero slope of the cor- 
relators at their onsets but revealing no additional effect. 
By this means, as in the shot-noise spectrum, the depen- 
dence on the bias voltage reveals again the onsets of the 
four correlators. Since we have chosen positive hQ^ the 
auto-terminal contributions are non-zero over the whole 
bias range. As before, cross-terminal ones vanish if no 
energy is provided by the voltage source: Cl^r 7^ if 
eV > —hft and Cr^l ^ if eV < hft. As it is clear from 
the bottom plot of Fig. 13, the oscillations of the com- 
ponents are not in phase, thus adding up to complicated 
oscillations in Sll{^)' But, as mentioned in the begin- 
ning, we doubt this could be a measurable effect. The 
shot-noise and the auto-terminal correlators are symmet- 
ric in the voltage dependence if Vq = 0, whereas the cross 
terminal ones obey Ca^fdi^^V) = C(3^a{^^—V). Here 
the charge-neutrality point and the width of the region 
III5 are revealed as a minima in the slope of the correla- 
tor Cl^l{^) at eV = :^2hQ and in the change of sign in 
dCR^R{n)/dV at eV = 0. 



VII. CROSS-CORRELATION NOISE 



The explicit expressions of Eq.(27) for the cross- 
correlation current noise spectrum of Fig. 14 can be ex- 
tracted from the general expression Eq.(21) in the same 
way as we did when deriving Eq.(24). From Figs. 8 a) 
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and b) it is also clear that the spectrum of auto-terminal 
correlators are oscillating as a function of Q with larger 
amplitude than cross-terminal ones, since they show an 
alternating behavior between positive and negative inte- 
grands. Dependence on q in the relevant frequency range 
is weak, as shown in Fig. 8 a). Contrary, cross-terminal 
contributions as in Fig. 8 b) show features with an alter- 
nating sign along both variables, Q and q. Thus, integra- 
tion along y-momentum leads to averaging and therefore 
significantly smaller oscillation amplitudes occur. As dis- 
cussed for the excess noise of the auto-correlation noise 
spectral function, we find the oscillations have a fre- 
quency hQz = 2eF what is tantamount to a period 
T = 27r in the plots. Complex conjugation corresponds to 
time-reversed states. Again, as the product of scattering 
matrices of the integrands in Eq.(21) suggests, it is prob- 
ing transmission- and refiection amplitudes of electron- 
hole pairs separated by an energy quanta HQ. So, for 
cross-terminal noise not only the refiection but also the 
complex transmission amplitude is essential even with- 
out ac-bias voltages. Again it would be interesting to 
test if the resulting oscillations could be detected in the 
challenging task of a finite-frequency cross-correlations 
experiment. Analyzing the integrands reveals the sym- 
metry C^^^{Q) = C'^^^{Q) if fiL = —fJ^R as we show in 
Fig. 15. This symmetry is distorted by applying an off- 
set voltage Vq. The spectrum of the correlator (7£^^(r^) 
shows a shift of the maxima and minima of the oscilla- 
tions with respect to C^^^(l^) for finite Vq. This shift 
is due to the fact that the distance between neighboring 
maxima of the integrand is not constant when varying 
hft at given q-mode (see the bending of the maxima to- 
wards higher frequencies for larger q in the integrands, 
e.g. Fig. 6). Derivatives of the correlators C^^o{^) with 
respect to voltage show a sequence of pairs of different 
maxima. This observation is traced down to the same 
origin as above, and so the appearance of peculiar os- 
cillations in the summed up cross-correlation shot-noise 
Slr{^) is explained. At hQ = current conservation and 
the unitarity of the s-matrix require Slr{^) = — ^ll(^)- 
Therefore the correlator described by Eq. 27d is negative. 



VIII. FINITE-FREQUENCY NOISE AT 
AC-BIAS 



By applying an ac-bias voltage at the leads one can 
inject charge-carriers at positive and negative energies of 
the Dirac cone without applying a dc- voltage. Analogous 
to the minimal conductivity, in the non-driven case going 
along with a maximal Fano factor, the shot-noise at zero 
frequency but finite ac-bias Saa{^ = 0;co') mirrors the 
behavior of the conductivity in Fig. 3. The noise spec- 
trum Fig. 16 for the driven setup (a = 1) is similar to 
the one without driving but with additional steps in the 
derivatives. For arbitrary ac-bias these steps can appear 
at frequencies hQ = {jUa — /i^) ± nhw due to the onset of 
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FIG. 15: (color online) Real parts of current-current cross- 
correlations in units of 27r/i/e^, as function of symmetrically 
applied dc-bias voltage and for fixed frequency. Jumps in the 
derivatives at ±QL/vf are due to the onsets of the cross- 
terminal contributions. 
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FIG. 16: (color online) Derivatives of current- current correla- 
tions real parts in units of 27Th/e^ with respect to frequency 
as function of frequency. We have chosen a symmetrically 
applied dc-bias voltage with additional harmonic ac-driving 
(ujL/vf = 4, a = 0.5 and a = 1) in lead L. Dashed vertical 
lines mark the step-positions, coloring specifies the correlator 
which shows the step at corresponding Q.L/vf- 



higher-order PAT events. Since we set a = 1 in Fig. 16, 
the correlator with states R ^ R shows no ac-induced 
steps in the derivative. But when \a\ ^ 1 all integrands 
(Fig. 9) are not given in terms of probabilities and can 
take negative values as mentioned in section V. For the 
shot-noise spectrum there are then two possible sources 
of contributions that could reduce noise: Either a corre- 
lators integrand or the the product of Besselfunctions is 
negative. When the driving voltage is applied symmetri- 
cally (a = 0) more PAT-induced steps in the derivatives 
of the noise spectrum are visible and finite contributions 
at negative Vt are possible for all correlators. As proposed 
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by Trauzettel et. al.,^-^ a time-dependent voltage could 
be used to induce interference between states in particle- 
and hole-like parts of the Dirac spectrum. This should 
correspond to Zitterbewegung like in relativistic quantum 
mechanics, but we are not aware of any unique feature 
caused by Zitterbewegung that can be distinguished from 
other oscillations, especially of Fabry-Perot nature. 



IX. CONCLUSIONS 

We have analyzed conductivity and non-symmetrized 
finite- frequency current-current correlations for a Fabry- 
Perot graphene structure. Oscillations on the intrinsic 
energy scale L/hvp are still present in the finite fre- 
quency noise. Emission spectra are diverging for large 
frequencies, whereas the absorption branch of the spec- 
trum has to vanish at h£L = —eV. As expected from 
the integrands, the current-noise also diverges for volt- 
ages \eV\ ^ hvp/L. Since the onset of the different 
noise contributions is defined by the four possible combi- 
nations of the chemical potentials, the noise built by all 
correlators consists of contributions oscillating with the 
same period but different phases. Although dominated 
by Cl^l(^) when correlating the currents in terminal 
L at large frequencies, this interplay is revealed in the 
spectra and voltage dependence of all correlators. Each 
contribution can show peculiar oscillations at low enough 
frequencies or voltages. In this regime features in the in- 
tegrands {q^ e)-dependence can have a prominent impact 
whereas they tend to be averaged out at large frequen- 
cies. Another aspect is the appearance of a special region 
showing phase jumps in the energy dependence of the in- 
tegrands when hTt < 2eV. This interplay of the Dirac 
spectrum and the Fabry-Perot physics^^'^^ can be probed 
purely by applying an appropriate combination of dc-bias 
and offset voltage Vq = F/2, thus connecting electron- 
and hole part of the Dirac spectrum symmetrically when 
eV = 2h£t. The way the scattering amplitudes are com- 
bined in this approach spoils the clear picture in terms of 



transmission- and refiection probabilities. Instead, in the 
dc-limit it gives rise to the interpretation of the L ^ L 
contribution in terms of jumps in the scattering-phase 
between time-reversed electron-hole states separated by 
the photon energy hQ. In the same way the complex cor- 
relators for cross-correlation noise or for the driven setup 
exhibit phase jumps and can not be written in terms 
of probabilities. Complex contributions of the scatter- 
ing matrices lead to large oscillations between positive 
and negative values of cross-correlation noise or in the 
derivatives with respect to frequency of the auto-terminal 
noise spectral function. These oscillations have a fre- 
quency of hQz = 2eF, what corresponds to a period of 
T = 27r in our plots. This frequency corresponds to the 
Zitterbewegung frequency as it is known for relativistic 
Dirac fermions. Again, strongly non-harmonic features 
can occur when the transition between different regimes 
is probed, especially when region III5 around the Dirac 
point comes into play. Additional ac-bias complicates the 
picture because combinations of q^ hVt^ mhuj define addi- 
tional phase jumps, onsets of the correlators and there- 
fore steps in the noise when higher-order PAT events oc- 
cur. Then the special role of the complex refiection and 
transmission amplitudes is essential for all possible cor- 
relators. 



Acknowledgments 

We would like to acknowledge the financial support by 
the DFG (Grant No. SFB 767) and thank B. Trauzettel 
for validating and sharing corrections to reference^^. 



Appendix A: Noise formulas 

The non-symmetrized noise spectrum under harmonic 
ac-driving without interactions is determined by 



5«,(^,^) = fe) 



de ^ Ji 



^6,lkm 



Jk I —^ I Ji 



m-\-k — l I :.- I ^rr 



eK< 



njjj J \ Rjjj J \ huj J \ huj 

Tr [A^5{a, e, e + hn)As^{p, e ^ hft ^ {m - l)hw, e + (m - l)hjO, )] f^{e - Ihu) (1 - fs{e ^ HQ - khu)) . 

I 



(Al) 



For the Fermi distribution function in lead 7 we use the 
shorthand /7(e) = l/(exp [(e — /i^)//cBT] + 1). In the 
limit of ksT — the distribution functions in the leads 
are given by Heaviside-Theta functions 0(/i^ — e) that 
define the integration intervals. Explicitly writing down 
the expression above for chosen Qf,/3 then leads to the 
four possible contributions to auto-correlation noise of 



Eq. (24) and to cross-correlation noise of Eq. (27) via 
summation over reservoir indices 7, 5 = L^R. 
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Appendix B: Boundary conditions 

As in^^ we confine the charge-carriers along the y- 
directions by infinite mass boundaries that diverges at 
the edges ^ = 0, W . This corresponds to the boundary 
conditions'^^ 



The solution is straightforward and determined by the 
four complex coefficients a(e), 6(e), t(e), r(e). The trans- 
mission t{e) and refiection r(e) amplitudes define the s- 
matrix and thus the current and noise of the scattering 
device at all voltages. 



^1 



^5 



y^Q 



y^O 



^1 



^, 



=W 



^W 



(Bl) 



Unlike the procedure for the Schrodinger equation, in 
graphene one only has to match the wave function itself 
and no constraint is given for the derivatives. Now ex- 
ploiting the boundary conditions along the x-direction, 
transmission- and refiection amplitudes are fully deter- 
mined by these constraints of the field operators at the 
Fermi levels. Therefor it is sufficient to match the wave 
functions at x = 0, L without ac-driving. A plain wave 
ansatz to solve the Dirac equation (1) for an electron 
incident from the left {x < 0) with energy e is given by 



^(x) = < 



%\ 



r(e) ^o^^_ if X <0 
«(e) ^S!V + Ke) ^S,- if 0<x<L (B2) 



if X > L 



where /i:(e),/^^(e) are the complex wave vectors in the 
reservoirs and and k{e) the one in the sandwiched 
graphene strip. Matching conditions at boundaries (con- 
tinuity at X = 0,L) combined with high doping in the 
reservoirs lead to the following set of coupled, complex 
equations: 



l + r(6) e 

1 - r(e) 



.(e)/2 



b{e)- 



^ia(e) 



Y^cos(a(e)) Y^cos(a(e)) 



V2 
tie) 



a{e) 



piaie)/2 



b{e)- 



-ia{e) 



i/cos(a(e)) i/cos(a(e)) 

5(e) 



(B3) 



(B4) 



gi/c(e)Lg-ia(e)/2 



^ = a(e) 

v2 ycos(a) 



g-i/e(e)LgiQ;(e) 



Y^cos(a) 



(B5) 



gi/c(e)Lgia(e)/2 



v2 ycos(a) 



6(e) 



g-i/c(e)Lg-iQ;(e) 



Y^cos(Qf(e)) 



(B6) 



t(e) = 
r(e) = 

a(e) = 

5(e) = 



1 



cos(fc(e)L) — i sec(a(e)) sin(fc(e)L) 

1 
— cot(fc(e)L) cot(a(e)) + i csc(a(e)) 

\/2cos(^) Vcos(a(e)) 

1 + e2ifc('=)^(-l + cos(a(e))) + cos(a(e)) 

eifc(Oi^cos(a(e))sin (^) 

V2(i cos(fc(e)L) cos(a(e)) + sin(fc(e)L)) 



(B7) 
(B8) 

(B9) 

(BIO) 
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